Self-Excitation and Feedback Cooling of an Isolated Proton 
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The first one-proton self-excited oscillator (SEO) and one-proton feedback cooling are demon- 
strated. In a Penning trap with a large magnetic gradient, the SEO frequency is resolved to the 
high precision needed to detect a one-proton spin flip. This is after undamped magnetron motion is 
sideband-cooled to a 14 mK theoretical limit, and despite random frequency shifts (typically larger 
than those from a spin flip) that take place every time sideband cooling is applied. The observations 
open a possible path towards a million-fold improved comparison of the p and p magnetic moments. 



The demonstration of a one-electron self-excited os- 
cillator (SEO) p] and one-electron feedback cooling [2] 
eventually led to greatly improved measurements of the 
electron magnetic moment and the fine structure con- 
stant [3] . An electron spin flip caused the SEO oscillation 
frequency to shift observably because of a magnetic gra- 
dient added to the Penning trap that held the electron. 

This Letter reports the first one-proton SEO and the 
first one-proton feedback cooling. The SEO frequency is 
resolved at the very precise level needed to observe a spin 
flip of a single p or p - opening a possible way to compare 
the p and p magnetic moments [3 H] a million times 
more precisely than current comparisons. To compensate 
for the p and p moments being 650 times smaller than 
that of an electron, these observations took place in a 
smaller trap with a much larger gradient than sufficed 
for observing spin flips of free [3] and bound [5] electrons. 
Feedback cooling promises to narrow the spin resonance 
linewidth and increase the transition rate. 

Effective sideband cooling [5J [7J is essential for these 
measurements, to center the p or p so that electric and 
magnetic anharmonicities do not mask the tiny frequency 
shifts that would signal spin flips. Sideband cooling is 
widely utilized not only in precision measurements of 
trapped particles, but also for the manipulation of qubits 
[5], the most precise clocks [5], the cooling of neutral 
atoms in lattices [TU] , and the cooling of mechanical res- 
onators [IT]. The one-proton SEO and feedback cool- 
ing make it possible to directly investigate the outcome 
of sideband cooling to the low-temperature, theoretical 
limit [7J Q~2] , for a particle that has no internal degrees of 
freedom available for cooling. The observed distribution 
of proton magnetron states is at the 14 mK theoreti- 
cal limit. Previous such investigations with an electron 
[7J [T3] reached only an energy 400 times higher than the 
desired low temperature limit, and investigations with 
an ion probed temperatures orders of magnitude higher 
[21 [TS] . A complication due to the strong magnetic gra- 
dient required is that an application of sideband cooling 
randomly shifts the SEO frequency by an amount larger 
than the shift from a spin flip on average. 

A single proton is suspended in a vertical B = 5.68 
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Tcsla field at the center of a cylindrically symmetric trap 
(Fig. [T]) - stacked rings with a 3 mm inner radius. The 
electrodes and surrounding vacuum container are cooled 
to 4.2 K by a thermal connection to liquid helium. Cry- 
opumping of the closed system made the vacuum better 
than 5 x 10 -17 Torr in a similar system [16], so colli- 
sions are not important. Appropriate potentials applied 
to copper electrodes (with an evaporated gold layer) in an 
open-access geometry [TTJ make a very good electrostatic 
quadrupole near the trap center, while also maintaining 
an open access to the trap interior from either end. 
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FIG. 1. (a) Penning trap electrodes and radiofrequency 
schematic for feedback cooling and self-excitation of the pro- 
ton axial motion. The feedback has gain G and a phase shifted 
by <f>. (b) Energy levels and transitions (arrows) involved in 
axial sideband cooling of proton magnetron motion. 

The proton's circular cyclotron motion is perpendic- 
ular to B, with a frequency uj + /(2tt) = 79.5 MHz that 
is slightly modified by the electrostatic potential. The 
proton also oscillates parallel to B at uj z /(2tt) = 553 
kHz. Because the potential is not a perfect quadrupole, 
this frequency depends slightly upon oscillation ampli- 
tude, A, with uj z {A) ss uj z . The proton's third motion is 
a circular magnetron motion, also perpendicular to B, at 
the much lower frequency, w_/(2/t) = 1.9 kHz. 

To couple the proton spin moment and the magnetron 
state that is the outcome of sideband cooling to the mea- 
surable uj z (A), the trap's central ring electrode is made 
of saturated iron (unlike the copper endcap and compen- 
sation electrodes above and below) . The extremely large 
magnetic bottle gradient, 



AB = /3 2 [(z 2 -p 2 /2)z-zpp} 



(1) 



with /3 2 = 7.8 x 10 4 T/m 2 , is 51 and 8 times larger than 
what was used to measure free [3j and bound [51 [T3] elec- 
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tron magnetic moments. The bottle reduces the field 
within the trap by 0.47 T (8 %). 

The axial frequency uj z (A) depends primarily on the 
strength of the z 2 term in the electrostatic quadrupolc. 
A magnetic moment /iz (from circular cyclotron or mag- 
netron motions, or from spin) adds a term going as jiz 2 
to the trapping potential, shifting uj z (A) by 
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The magnetic moments from cyclotron and magnetron 
motion go as n and £. The 553 kHz axial frequency shifts 
by 21 mHz per cyclotron quantum (200 Hz per fim for our 
typical cyclotron radius), and by 0.5 /iHz per magnetron 
quantum (40 mHz per /im at our typical magnetron ra- 
dius). A proton spin flip will cause a 60 mHz shift. 

The frequency of the current induced to flow through 
R in Fig. [T] by the axial oscillation is measured to deter- 
mine uj z (A). The voltage across R = 25 Mfi is Fourier 
transformed after an amplifier that uses a high electron 
mobility transistor (HEMT) with good thermal connec- 
tion to 4.2 K. The I 2 R loss for the induced current going 
through R gives an axial damping time 7J 1 = 60 ms. 

R represents losses in an LC tuned circuit resonant at 
U! z . The losses are minimized to maximize R, the ob- 
served signal and the damping rate. Varactors tune the 
circuit and its matching to the HEMT. A superconduct- 
ing inductor with L = 2.5 mH cancels the reactance of 
the trap capacitance, leaving R = Qui z L. The circuit's 
quality factor is tuned to Q = 3000 with the trap elec- 
trodes attached. (With the amplifier connected only to 
the endcap, Q = 5000. With a capacitor substituted for 
the trap electrodes, Q > 20, 000.) The proton's cyclotron 
and magnetron motions in this trap are not damped. 

A nearly identical "precision trap" is just below the 
trap in Fig. [I] A detection circuit resonant at uj + / {2tt) — 
86.5 MHz, attached across halves of a copper ring elec- 
trode, damps this motion in 77 = 10 min. (This could 
be three times faster with better amplifier tuning.) An 
axial amplifier detects and damps the axial motion. 

A single proton is isolated in the second trap using a 
relativistic method we developed earlier with antiprotons 
|18j . An H atom is ionized in the trap by an e~ beam 
from a sharp field emission point. Strong driving forces 
applied at the axial frequencies of unwanted ions keep 
them from loading. A strong pulse of cyclotron drive 
produces one-proton cyclotron resonances that differ in 
frequency because of differing cyclotron energies and rel- 
ativistic mass shifts of w+. The trap potential is tem- 
porarily reduced until the signal from only one proton 
remains. The cyclotron energy of the remaining proton 
damps until its radius is less than the 0.5 /im average 
for a 4.2 K distribution. After magnetron cooling, the 
proton is transferred into the trap of Fig. [l] by adjusting 
applied potentials to make an axial potential well that 
moves adiabatically from the lower to the upper trap. 

The proton axial oscillation whose frequency is to be 



measured satisfies the equation of motion, 

z + ~/ z z + [u z (A)+Au z } 2 z = F d (t)/m. (3) 

A driving force Fd(t) is added to the restoring force (from 
the electrostatic quadrupole and the magnetic bottle), 
and to the damping force —m^f z z (from the loss in R). 

With no feedback, F d is the Johnson noise from the 
resistor that is then amplified and detected. The proton's 
axial oscillation shorts this noise 19J, making a dip in the 
noise power spectrum (Fig. ^p) whose half width is j z . 
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FIG. 2. SEO peak (a) and noise dip (b) for 160 s of averaging, 
(c) Frequency resolution achieved with a single average of an 
SEO peak (black x) and noise dip (red x), with the standard 
deviation (black points) and Allan deviation (black triangles) 
of averaged SEO measurements, (d) Drift of 256 s averages 
over sixteen night time hours. 

The axial frequency is determined to higher precision 
using the better signal-to-noise and narrower signal width 
of a one-proton SEO (Fig. [2^,). The one-particle SEO, re- 
alized previously only with an electron pQ, is realized by 
adjusting the amplitude and phase of the amplified in- 
duced signal and feeding this back to the other side of 
the trap as a driving force on the proton. The feedback 
produces a force Fd(t) ~ mGj z z with feedback gain G. 
Self-excitation occurs, in principle, when the feedback 
cancels the damping at G = 1. Noise causes amplitude 
diffusion and energy growth, however, and G slightly dif- 
ferent from unity will either decrease or increase A expo- 
nentially. A stable and useful SEO thus requires limiting 
the amplitude to some value A Q . Here a digital signal 
processor (DSP) chip Fourier transforms the signal to 
determine A, and makes G = 1 + a(A — A Q ) [I]. 

An axial oscillation z(t) — Acos(u>t) generates a feed- 
back force Fd(t) = —uAGm^ z sin(wt + </>), when a phase 
shift <f> is introduced (Fig.[T]). Inserting in Eq. [3] yields 

Gcos(0) = l (4) 

Gw lz sm^)=Lo 2 -\lo z {A)] 2 (5) 

For ("f z /uj z ) tan(c/)) << 1, the SEO thus depends on tj> as 



w(A,4>) « uj z (A) + ^tan(4>) 



(6) 
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With positive feedback, and a feedback phase adjusted 
to optimize the signal (Fig. [3^,) , the measured SEO fre- 
quency as a function of feedback phase fits well to Eq. [6] 
(Fig. p5p). The scatter in repeated frequency measure- 
ments(Fig. [3J;) is reduced when the trapping potential is 
tuned to make the best possible electrostatic quadrupole. 
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FIG. 3. (a) SEO signal strength vs. feedback phase, (b) 
Measured axial frequency vs. feedback phase (points) fit to 
the expected Eq. [Fj] (c) Tuning for optimal SEO stability by 
adjusting trap anharmonicity. 

Sideband cooling, a method to radially center the pro- 
ton, is especially important given that the magnetic field 
changes significantly as a function of radial position. A 
sideband cooling drive at uj z + w_, applied across the 
halves of a compensation electrode (Fig. [TJ to produce 
an zy potential gradient [7J, makes the transitions be- 
tween axial states k and magnetron states I indicated by 
arrows in Fig. [TJd. The probability Pj~.i satisfies 



= -P k ,t 



-Pui+i T- + ^fc+ii-i r+- (7) 



This steady-state rate equation has rates that depend on 
axial and magnetron raising and lowering operators [7J, 

T+ ~ \ <k + l,t-l\aia-\k,£>\ 2 ~ (k + l)£ (8) 
T_ - |<fc- 1,1 + l\a z al\k,£>\ 2 ~ k(£+l). (9) 

The axial distribution remains a Boltzmann distribution 
due to its coupling to the detection resistor, a reservoir at 
temperature T z , so that Pk t i = pt exp [— ktlD z /(kBT z )]. 
The solution to Eq. [7] is a magnetron distribution, 



pi ~ exp [-ihuj- / (k B T m )], 



(10) 



where the effective magnetron temperature is T m — 
T z u}-/uj z . The theoretical cooling limit [71 112). 

k B T m = (-E mag ) = — (E z ) = — k B T Zl (11) 

0J Z UJ Z 

comes from evaluating the average of the magnetron en- 
ergy, (E mag ) = 52 e PiE e , with E t = -(£ + l/2)fiw_. 
Decreasing £ increases the potential energy on a radial 
hill, while decreasing the much smaller kinetic energy, so 
Ei is negative. Ref. |7j extends the theoretical argument 
and limit to off-resonant sideband cooling drives. 

The outcome of sideband cooling in the trap of Fig. [I]} 
is investigated with a three-step sequence. First, the ax- 
ial energy is either left in equilibrium with the detection 
resistor or modified using feedback. Second, a sideband 
cooling drive at uj z + u;_ is applied and then turned off. 



Third, the SEO is started and the axial frequency mea- 
sured. Each application of sideband cooling produces a 
measurably different magnetron state and uj z (A), so the 
sequence is repeated to make histograms of measured ax- 
ial frequencies. 

The gray histogram and Gaussian fit in Fig. [4^, show 
the scatter for repeated measurements of u> z (A) taken 
with no sideband cooling drive (i.e. no change in mag- 
netron radius) and no feedback (i.e. no change in T z ). 
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FIG. 4. (a) Histograms of magnetron states after no sideband 
cooling (gray), and produced by sideband cooling using feed- 
back cooling (blue), no feedback (green), and with feedback 
heating (red) . Solid curves are convolutions of the gray Gaus- 
sian resolution function and Boltzmann distributions at the 
specified T z . (b) The magnetron radius increase from a side- 
band drive at uj z — ui~ is fit to an exponential and extrapolated 
back to an initial magnetron radius. 

Sideband cooling with no feedback broadens the gray 
into the green histogram (Fig. [4^,) . A convolution (green 
curve) of Eq. 
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with T m = 30 mK (corresponding to 
T z = 8 ± 2 K) and the gray Gaussian resolution function 
fits the measured histogram when Eq.[2]is used to convert 
magnetron energy to axial frequency shift. The axial 
temperature is reasonably higher than the T z = 5.2 K we 
realized with one electron in a 1.6 K apparatus [2]. 

Feedback changes the measured T z as predicted, from 
T z0 at G=0 to T Z (G) = (1 - G)T z0 (Fig. [5Jd), increas- 
ing our confidence in this new way to measure the axial 
temperature. The corresponding damping widths also 
change just as predicted, from T z0 to T Z (G) = (1— G)T z0 
(Fig. [5]i). The ratios in Fig. [BJ: are constant, consistent 
with the fluctuation-dissipation theorem. 

Feedback cooling of the axial motion to T z = 4 K nar- 
rows the distribution of magnetron states (blue histogram 
in Fig. [4^,) such that the effective magnetron tempera- 
ture is T m — 14 mK. Feedback cooling from T z — 8 to 
4 K seems plausibly higher than the cooling from 5.2 to 
0.85 K [2] achieved with one electron in a 1.6 K appara- 
tus. Feedback heating of the axial motion to T z = 20 K 
broadens the distribution (red histogram in Fig. |4^i) . 

A check on the magnetron orbit size produced by side- 
band cooling comes from expanding the orbit size expo- 
nentially with a sideband heating drive at lu z — w_ [71ll3j. 
Each trial (e.g. Fig. [4J5) is extrapolated to determine the 
radius at the start of the heating. Averaging the ini- 
tial radii from 200 trials gives 11 ± 2 fim, consistent with 
T z = 8±2 K from Fig. [4^,, and hence with the theoretical 
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FIG. 5. Measured damping widths (a), temperatures (b), and 
their ratios (c) as a function of the feedback gain, G. 

cooling limit (Eq. [T7|). We do not understand the earlier 
electron observations [7J [T3] but note that progress has 
been made in the detection electronics that sets T z . 

Fig. [2J; shows that for the first time in a strong mag- 
netic gradient, the resolution achieved in measuring u z is 
comparable to the 60 mHz resolution needed to observe 
a proton spin flip (for a 16 s averaging time). The resolu- 
tion is much better for a single SEO measurement (black 
x) than for a single dip measurement (red x). Repeated 
SEO frequency measurements have a standard deviation 
(black points) and an Allan deviation (black triangles) 
that is larger than the precision for a single measure- 
ment. Fluctuations in the trapping potential, mechanical 



vibrations, temperature variations, and fluctuating con- 
tact potentials are being investigated as possible sources 
of the scatter. 

Sideband cooling is required to minimize the scatter 
and achieve the 60 mHz resolution. However, Fig. [4^ 
demonstrates that an application of sideband cooling ran- 
domly selects a new magnetron radius and u z that is 
shifted by much more than this resolution on average. A 
solution starts with an initial sideband cooling period, 
after which uj z is measured. An attempt to make a spin 
flip can be then made for several minutes, during which 
time unwanted magnetron heating shifts u> z typically by 
about 0.3 Hz/hr (Fig.[2]i). Measuring lo z will thus reveal 
a shift larger or smaller than 60 mHz depending upon 
whether the spin has or has not flipped. The process can 
then be repeated, with cooling selecting a new lu z . 

In conclusion, a one-proton self-excited oscillator and 
one-proton feedback cooling are realized for the first time. 
A very strong magnetic gradient is added to the Pen- 
ning trap in which the proton is suspended to make it 
possible to observe sideband cooling distributions and to 
investigate the possibility of observing spin flips. Side- 
band cooling of the undamped proton magnetron motion 
to 14 mK is demonstrated, even though every applica- 
tion of sideband cooling shifts the monitored SEO os- 
cillation frequency more than would a spin flip. As an 
application, the SEO oscillation frequency is resolved at 
the high precision needed to observe a spin flip of a single 
p or p, opening a possible new path towards comparing 
the p and p magnetic moments at a precision higher than 
current comparisons by six orders of magnitude or more. 
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